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EXACT ANALYTICAL SOLUTION TO A TRANSIENT CONJUGATE 
HEAT-TRANSFER PROBLEM 
by James Sucec* 

Lewis Research Center 

SUMMARY 

An exact analytical solution is found for laminar, constant- property, slug flow over 
a thin plate which is also convectively cooled from below. The solution is found by 
means of two successive Laplace transformations when a transient is initiated by a step 
change in the fluid inlet temperature. Two cases are considered: (1) the case where 
the fluid and the plate are originally at the coolant temperature before the step change 
in inlet temperature and (2) the case where the fluid and the plate are at the eventual 
steady state of the first case when another step change in inlet temperature occurs. 
Results are given for the temperature within the moving fluid as a function of space 
coordinates and time, for wall surface temperature, and for surface heat flux. The re- 
sults of the exact solutions are compared with the quasi- steady values, and a criterion 
for the validity of the quasi- steady solution is evolved. Also the effect of the fluid- to- 
plate coupling is investigated. 

One notable difference between this type of transient (a step change in the fluid inlet 
temperature) and the transient induced by a step change in wall temperature is the ab- 
sence of infinite and very large heat fluxes, which are characteristic of the transient 
induced by a step in the wall temperature. 


INTRODUCTION 

This report analytically predicts the transient surface heat flux and temperature, 
on a cooled turbine blade or vane, caused by an abrupt change in the inlet temperature 
of the gas stream . 


* Associate Professor of Mechanical Engineering, University of Maine, Orono, 
Maine; Summer Faculty Fellow at the Lewis Research Center in 1972. 



In the design of blades and vanes for gas turbine engines, situations often arise 
where the fluid flowing over a solid boundary undergoes a transient, thus inducing time- 
dependent surface temperatures and heat flux. More often than not, a quasi-steady 
analysis is applied to such a problem; that is, it is assumed that the steady-state rela- 
tions for the heat flux and the surface coefficient of heat transfer are approximately 
valid at each instant of time so long as the instantaneous values of temperature and 
velocity are inserted into these relations. That the quasi-steady approximation can in- 
duce appreciable error has been demonstrated in many publications, among them refer- 
ences 1 to 6. Sparrow (ref. 1) solves the problem of a step jump in the surface temper- 
ature at a two-dimensional planar stagnation point for laminar, constant-property flow 
when the velocity field is steady. By inserting Kantorovich velocity and temperature 
profiles into the integral form of the energy equation, he obtained approximate solutions 
for the transient surface heat flux and surface coefficients of heat transfer. In particu- 
lar, he found that for small times, the process was that of pure conduction, with the 
result that the heat flux at the surface is inversely proportional to the square root of the 
time from the step change in temperature. Hence, very large (relative to steady- state 
values) surface heat fluxes are present during the initial portion of this type of transient. 
The linearity of the constant- property energy equation allows the response to the step 
change in wall temperature to be generalized to a wall temperature which varies arbi- 
trarily with time. This he accomplished by way of the superposition integral. Siegel 
and Sparrow (ref. 2) solve for the transient surface heat flux due to a step change in 
wall temperature for the laminar thermal entrance region of flat ducts. Again, Kantoro- 
vich approximating sequences are used for the velocity profile (steady) and for the un- 
steady temperature profile. This time, however, insertion of these profiles into the 
integral form of the energy equation yields a partial differential equation for the thermal 
boundary layer thickness as a function of axial distance and time. They use the method 
of characteristics to arrive at the solution, which is then generalized by superposition 
to handle wall temperature variations that are arbitrary functions of time but are spa- 
tially uniform. The step change in surface flux and its generalization is also considered. 
The problem of transient heat transfer for laminar slug flow in ducts yields to the tech- 
nique of separation of variables in an exact analytical solution by Siegel in reference 3. 
Cess (ref. 4) solves for the transient surface heat flux in laminar, constant-property, 
steady- velocity flow over a flat plate which has its surface temperature changed ab- 
ruptly. He finds a series solution for small time in the physical plane, and one for long 
time in a Laplace transform plane. These limiting. solutions are then joined, in an ap- 
proximate manner, in the Laplace transform plane, and the result is inverted to give the 
surface flux over the entire time domain. Goodman (ref. 5) uses a Lighthill-type veloc- 
ity profile in the integral energy equation to yield an equation which .can be solved by 
characteristics for arbitrary free- stream velocity variation with position along the sur- 
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face of interest when the wall temperature abruptly changes. In reference 6, Chao and 
Jeng deal with unsteady heat transfer due to either wall temperature variation or wall 
flux variation with time at either a two-dimensional planar or an axisymmetric stagna- 
tion point. Sparrow and Siegel (ref. 7) consider unsteady turbulent heat transfer in 
tubes due to a step change in the wall temperature at time zero. Chambr§ (ref. 8) uses 
a double Laplace transformation to arrive at an exact analytical solution for the tran- 
sient heat transfer in laminar slug flow of a fluid over a plate insulated on the bottom 
but containing an exponential generation term. Axial conduction in the plate is neglected 
and the temperature distribution in the transverse direction is lumped. Thus, the plate 
has no internal resistance to energy flow in that direction, but it does have finite ther- 
mal storage capacity. This finite thermal storage capacity couples the transient tem- 
perature distribution in the fluid to that of the plate and, hence, it becomes a transient 
conjugate problem (see Perelman, ref. 9). Transient local surface temperatures of 
the plate are presented for the case where both the plate and the fluid are at zero non- 
dimensional temperature initially, then the exponential generation term is turned on 
while the fluid approaching the leading edge of the plate stays at a nondimensional tem- 
perature of zero. Note is taken of the plane "wave” phenomenon peculiar to slug flow, 
namely that the "front” of fluid which is located at the leading edge, at the instant when 
the generation in the plate is turned on, separates the thermal fluid field into two vastly 
different regions. Ahead of this front, the process is one of conduction in the fluid with 
no mixing; whereas, behind the front there is mixing due to fluid being processed over 
the leading edge. Siegel (ref. 3) also discusses this effect in some detail. In refer- 
ence 10, Soliman and Chambre use double Laplace transformations to find an exact 
analytical solution for the case of laminar, constant- property flow with a steady veloc- 
ity distribution which is a linear function of distance away from a plate which undergoes 
either a step change in the surface temperature or in the heat flux. 

In reference 11, Inouye and Yoshinaga use an approximate integral method to solve 
two problems of unsteady heat transfer at a stagnation point. In their first problem, the 
wall temperature is constant and the transient is caused by a step change in the fluid 
temperature. They indicate that there is no published work on this type of problem. 
However, it seems to this author that, near the stagnation point, there should not be any 
difference between the solution for a step change in wall temperature and their solution 
for a step change in the fluid temperature. This is so, because when the fluid tempera- 
ture abruptly changes, it changes for all the fluid. There is no waiting for free- stream 
fluid which has had its temperature changed to be carried to any point on the body. 
Hence, the physical situation looks just like a fluid with constant temperature in contact 
with a wall, at a stagnation point, which has just had its temperature abruptly raised. 
Inspection of the relevant partial differential equations and the boundary and initial con- 
ditions also bears out this conclusion. As a final check, this author compared, for a 
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Prandtl number of unity, a few selected predictions for the surface heat flux from the 
paper of Inouye and Yoshinaga (ref. 11) with those of an earlier paper by Chao and Jeng 
(ref. 6) and found excellent agreement between them. In the second stagnation- point 
problem worked in reference 11, the wall has its rear surface insulated and is initially 
at the temperature of the fluid when suddenly the fluid temperature is changed. Axial 
conduction was neglected in the plate and the plate temperature was lumped in the trans- 
verse direction, but the plate had finite thermal storage capacity. Hence, this is a 
transient conjugate problem, as was the problem in reference 8. Approximate results 
for the surface heat flux and wall temperature are given as a function of time for vari- 
ous values of a transient coupling parameter. Lyman (ref. 12) solves for the wall heat 
flux and wall temperature at a stagnation point when the free stream suddenly exper- 
iences a step change in its temperature. His is a conjugate problem which he solves in 
an approximate fashion using Kantorovich profiles in the integral form of the energy 
equation for the fluid while retaining the partial differential equation for the solid wall. 
For the case of a semi-infinite wall an approximate solution for all values of time is 
found. Short- and long-time approximate solutions are found for a finite thickness plate 
insulated on its lower surface. Unlike the approach of reference 11, Lyman applies the 
step change in free- stream temperature to the outer edge of a thermal boundary layer 
and calculates the transit time needed for the wall to feel the effects of this change. 

Thus, for times less than the transit time, the wall flux is zero. The flux is never in- 
finite even at the time when the wall first learns of the free- stream temperature change. 
By that time a smooth continuous temperature profile has developed in the thermal 
boundary layer, thus negating the possibility of an infinite temperature gradient at the 
wall. 

Except for the conjugate problems of references 8, 11, and 12 and the present re- 
port, there does not seem to be any literature on transient conjugate problems. How- 
ever, some steady-state conjugate problems in forced- convection heat transfer can be 
found in references 9 and 13 to 17. 

The instigation for the present work was the transient in the fluid, flowing over gas 
turbine blades or vanes, induced by a change in the fluid inlet temperature due to startup 
or due to a change in the power level of an already operating engine. One wishes to pre- 
dict the transient surface flux and the surface temperature variation. This problem, 
unlike the other transient heat-transfer problems already solved in the references, is 
unique in the sense that the reason for the transient is not some prescribed wall tem- 
perature variation with time. Rather the transient occurs because of the change in the 
fluid temperature at the leading edge. In addition, a blade or vane does not usually have 
some prescribed temperature but rather has prescribed internal coolant conditions. 

And hence, the solid over which the fluid flows has its temperature distribution coupled 
to that of the fluid - a conjugate problem. 
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ANALYSIS 


A cooled gas turbine blade or vane will be idealized as a thin flat plate of thickness 
b which has its lower surface exposed to coolant at a constant temperature T c and 
which has a constant surface coefficient of heat transfer h c between the coolant and the 
lower surface. The hot combustor gases at temperature and velocity u^ flow 
over the top of the plate as depicted in figure 1. (All symbols are defined in appendix A. ) 


Gases at and n* 


y 

x 

v~ • — ~ 

Coolant at T c and h c 

Figure 1. - Schematic diagram of physical situation. 

The analysis is for laminar, constant-property, two-dimensional planar, boundary- 
layer-type flow without appreciable viscous dissipation. 



Solution for Step Change in Fluid Inlet Temperature When 
Fluid and Plate Are at Coolant Temperature 

At time t = 0 both the plate and the fluid flowing over the plate are at the coolant 
temperature T c when suddenly the fluid for x == 0 has its temperature changed to a 
value T m different from T c< (The problem of a step change in the inlet temperature 
after the steady state has been reached, which occurs when going to a different power 
level or in shutdown of the engine, can be easily handled once the solution to the case 
just described is available. ) The surface- heat- flux, plate-surface-temperature, and 
temperature distributions within the moving fluid are required for all values of time 
t > 0. In order to investigate the general character and trends of such a solution, some 
further idealizations are made which permit an exact analytical solution. These ideali- 
zations are a steady slug flow velocity field (which corresponds to vanishingly small 
Prandtl numbers), neglect of axial conduction in the plate (which can sometimes be 
justified physically on the basis that cooled blades and vanes are relatively thin walled 
and constructed of relatively low- thermal- conductivity material), and lumping of the 
temperature distribution in the y-direction in the plate (which can often be justified if 
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the plate Biot number, based upon b and the larger of the coolant- side surface coeffi- 
cient and the gas-side surface coefficient, is relatively small). With 0 = T - T^, the 
partial differential form of the thermal energy equation for the fluid is (see, e. g. , 
ref. 18) 



30 

30 

3 *6 

(1) 



+ u oo — = 

Q! f 

1 9 


at 

3x 


The boundary conditions are 





t = 0 x > 

0 

y > 0 

0 = *c = T c - T - 

(2) 

x = 0 


t > 0 y 

> 0 0 = 0 

(3) 


y — 00 t>0 x>0 0 is finite 


(4) 


Next make an energy balance on the plate material, neglecting axial conduction and ther- 
mal resistance in the y-direction. Note that the plate temperature at any x and t must 
equal the fluid temperature at the same value of x and t. Thus, the partial differential 
equation for the plate temperature becomes a boundary condition for equation (1), 
namely, 


y = 0 t>0 x > 0 M =-£. (0 - 0 ) + r-22 (5) 

3y kj c 9t 

where r = p C_ 0 b/k f and is a measure of the thermal storage capacity, per unit length, 
of the plate material. 

The solution of equation (1), subject to the side conditions (2), (3), (4), and (5), was 
approached by way of two successive Laplace transformations. First, the Laplace 
transformations of the equation and the appropriate boundary conditions with respect to 
time were taken. The transformed temperature function in the first transformed plane is 

0e' pt dt (6) 

Performing the indicated transformation yields 


e = j?[e] = r 
t-P Jo 


6 



(7) 


p?. 9 „ + u „M = af ^ 


0X 


3y 


x = 0 y>0 0 = 0 


y_ oo x>0 0 is finite 


( 8 ) 


y = 0 x > 0 — = — - — V r(P 0 - O 

3y k f V P/ 

Next, a Laplace transformation with respect to x is taken for equation (7) and the last 
two equations of the set (8). The transformed temperature function in the second trans- 
formed plane is defined as 


0 = j?[0] = 

x— s 



0e" sx dx 


(9) 


This operation yields the following ordinary differential equation and its associated 
boundary conditions: 


2 a A 0O s + p\ = 


d 0 l u o o 

dy 2 


0 = - 


( 10 ) 


at* 


set* 


y — oo 0 is finite 


y = 0 ^ =( ^ + rp^0 - 9 C (± + Hc 


dy \ k. 


f 


k f sp) 


( 11 ) 

( 12 ) 


Solving equation (10) and subjecting its solution to condition (11) and to condition (12), 
which is actually the conjugation condition, yield the solution for 0 as 


h r \ -y s+p)/a f 
u„ 0„( — + rp I e 

k f / 

0 = — —L - + 


e , 


p(u„s + p) 


Ko S + P h c 
+ — + rp ; 


s(us + p) 


(13) 


at* 


f 
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Since one is primarily interested in the surface-temperature distribution and the 
surface- heat- flux distribution, he can, in the spirit of Lighthill (ref. 19) invert equa- 
tion (13) with y = 0 and invert dl/dy at y = 0, rather than the entire temperature 
distribution. Since the successive inversion with respect to x and t of equation (13) 
is not all that formidable, the author has inverted equation (13) and then computed the 
wall temperature and surface heat flux by operations in the x,y,t- plane. As an internal 
check the transformed surface- heat- flux and surface-temperature distribution in the 
s, y, p- plane was also inverted and exact agreement was noted. All the inversions were 
performed by using tables of transforms (specifically, refs. 20 and 21) in conjunction 
with the substitution, translation, and convolution theorems of the operational math- 
ematics. Details are given in appendix B. 

When the nondimensional variables 

T = 

X 


2 





and 


e 



are introduced, the exact analytical solution for the nondimensional temperature excess 
ratio in the moving fluid becomes 


JL = 1 + u(r - 1) (-1 + erfc[Y] - e 2T?Y+7?2 {erf [e(r - 1) + r, + Y] -erf [77 + Y]}) 


(14) 


where u(r - 1) is the unit step function, that is, 
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u(t- 1) 


= 0 for r < 1 

= + 1 for t > 1 


The surface-temperature variation is found by setting Y = 0 in equation (14), which 
gives 


w 


= 1 + u(t - l)^-e^ {erf[e(r - 1 ) + 77 ] - erf[rj]}j 


( 15 ) 


The nondimensional surface heat flux is defined as follows: 


% = 



(16) 


The denominator in equation (16) is the surface heat flux that would obtain in steady, 
slug flow over an isothermal plate at temperature excess equal to e . 

It can be readily shown that 



With this and equation (14), the exact analytical solution for the nondimensional 
heat flux becomes 

Q w = e” ^ + 2 rje(T- 1 )] + |^ er f[ e ( T _ 1 ) + 77 ] - erf[rj]} for r 

= 0 for t < 1 

The exact analytical solution for the nondimensional temperature is seen to 
the side conditions of equation (1) and, as a final check on the solution technique, equa- 
tion (14) was substituted into the governing partial differential equation (1) and was found 
to satisfy it identically in the domain of definition of 0. Some limiting cases of equa- 
tions (15) and (17) have also been checked to ensure that they agree with physical reason- 
ing. If 77 — °° , interpreted as h £ — «> , the plate must be at the coolant temperature 
excess 0 Q , and the nondimensional flux should reduce to that appropriate to an 


surface 



J 


satisfy 


9 



isothermal flat plate, namely unity. By applying L'Hopital's rule to the second term of 
equations (15) and (17), it is readily seen that, as 

i - “ - e c % - i (is) 

The coupling parameter e - 0 corresponds to infinite thermal storage capacity of the 
plate. And since the plate is initially at 9 C , it would behave, with e - 0, as a heat 
reservoir at 8 c for all values of time. Hence, should be unity. Equations (15) 
and (17) properly reduce to this degenerate case also. 


Solution for Step Change in Fluid Inlet Temperature 
After Steady State Had Been Reached 

The eventual steady state reached in the previous case can be found by letting the 
nondime nsional time r — in equations (14), (15), and (17). This yields 


^5® = erfc[Y] - e 2TjY+r?2 erfcfa + Y] 


(19) 


= 1 - e 7 ? 2 erfc[Tj] 
0 c 

^ 77 2 

%,ss = V ¥lie erfc[7]J 


( 20 ) 


( 21 ) 


where Q^ y gs denotes the eventual steady- state value of and similarly for the 
functions B gg and 0 W> sg . 

Now the following problem is considered: the plate of figure 1 is in the steady state 
described by equations (19), (20), and (21), where the fluid inlet temperature is T^ , 
when suddenly the inlet temperature of the fluid is changed to a different value, Tq, 
which corresponds to an increase or a decrease of the power level of a gas turbine en- 
gine (a so-called acceleration or deceleration). Once again, we would like to predict the 
temperature distribution within the moving fluid and the surface- heat- flux and surface- 
temperature distributions in this type of transient. The governing partial differential 
equation is still equation (1) and side conditions (4) and (5) remain the same. Side con- 
ditions (2) and (3), however, take the following form: 
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t = 0 x > 0 y > 0 0=0, 


( 22 ) 


ss 

where 0„ o is the function of x and y given in equation (19), and 

SS 

x = 0 t>0 y > 0 0=T O -T oo (23) 

Because of the case previously solved, it is natural to seek a solution to this case by 
using the following decomposition: 

0(x, y, t) = (p (x, y, t) + 0 ss (x, y) (24) 

Inserting equation (24) into equations (1), (4), (5), (22), and (23) yields 

cp (x, y, t) = T q - T m + y(x, y, t) (25) 

where y(x, y, t) satisfies equations (1) to (5) if 0 c is replaced by T M - Tq. Thus, the 
solution for y(x, y, t) is equation (14) with 0 C replaced by - Tq and, of course, y 
in place of 0. Doing this and using equations (25), (24), and (19) yields the exact ana- 
lytical solution for the temperature distribution within the moving fluid for this type of 
transient as follows (details in appendix C): 

2 

A = erfc[Y] - e 2l?Y+77 erfcfo + Y] 

0 c 


T 0 " T °° 


T - T 

oo c 


u(t - 1)^1 + erfc[Y] - e 2r 7 Y+7 7 [erf[e(r - 1) + 77 + Y] - erf[Tj + Y]^ 


(26) 


Setting Y = 0 in equation (26) gives the surface temperature excess ratio as 


w 


2 /t_ - T \ 2 

1 - e^ erfc[7j] - u(t - 1)^— —j e 11 {erf[e(T - 1) + rj] - erfjrj]} (27) 


By operating on equation (26) and using equation (21), the transient surface heat flux 
can be put in the following convenient form : 
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A = 1 MYA.) for r> 1 

^W,SS V T oo " T C/V^W, SS/ 


% 

^w, s 


= 1 for t < 1 


For r > 1, ss is given by equation (21) while is given by equation (17). 


RESULTS AND DISCUSSION 

The first case considered is where both the fluid and the plate are initially at tem- 
perature T c and suddenly the temperature of the inlet fluid at x = 0 is changed to T M 
and held there. Plots are presented, for this case, of nondimensional surface heat flux 
Q w and nondimensional temperature excess ratio 0 W /# C in figures 2 and 3, respec- 
tively, as functions of a time parameter c(t - 1). The plots are presented for various 
values of rj ranging from rj = 0 (lower plate surface insulated) to 77 — °° (which corre- 
sponds to a constant plate temperature T ). The values of rj between these two ex- 
tremes were chosen not only to illustrate the trends of the flux and the temperature with 



0 .5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 


«(t- 1) 

Figure 2. - Heat -transfer response to step change in inlet temperature when fluid and plate are initially 
at coolant temperature. 
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Figure 3. - Wall temperature response to step change in inlet temperature when 
fluid and plate are initially at coolant temperature. 


77, but also to be in a reasonable range of possible application. A study of 77 shows, by 
its definition, that it is a measure of the coolant- side surface coefficient divided by the 
local gas- side surface coefficient. In the range of gas turbine application, this param- 
eter could vary from about 0 to unity, with 77 = 0 corresponding to a solid uncooled blade 
or vane and with tj = 1 corresponding to a relatively sophisticated internal cooling 
scheme such as impingement or impingement with crossflow. The choice of e(r - 1 ) as 
the abscissa rather than ct is dictated by the fact that because of the unit step function 
in the solution, is zero for r < 1 and only attains nonzero values for r > 1 . A 
similar statement holds for 0 W / 0 C - This is a result of the fact that the plate at position 
x does not realize the inlet temperature has changed until the fluid that was at the inlet 
at r = 0 reaches the position of interest. Inspection of figure 2 shows that as 77 de- 
creases there is a more pronounced variation of flux with time, as would be ex- 
pected based on the previous discussion of the limiting case 77 — 00 , in which there is 
really no transient in the plate . Study of figures 2 and 3 allows determination of the 
value that e(r - 1) must be greater than or equal to so that both the surface flux and the 
surface temperature excess ratio are within 5 percent of their steady-state values. 

This critical value of e(r - 1 ) is denoted as a* ; the values of a* are presented in the 
following table for various values of 77: 
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h c 

V =—'\ 
kf 1 

Ft 

1*00 

Value which a time parameter, e(r - 1), 
must be greater than for solution to 
be within 5 percent of steady state, 
a* 

0.25 

1.75 

.50 

1.375 

1.00 

.875 

OO 


0 


For t) = 0, this type of criterion does not hold since the steady- state values of and 
ejQ c are both zero. For this case, ij = 0, it may be noticed that for e(r - 1) 2. 2, 
both and 0 w /0 c are less than 0. 02. 

The influence of the coupling parameter e can be seen by reference to figure 4 
where is plotted versus r - 1 for different values of e at two different values of 
tj. One immediate observation is that the eventual steady- state is independent of e 
and dependent only on tj. That this should be so can be verified by recalling that neglect 
of axial conduction in the plate and the lumping of the plate temperature in the y-direction 
causes the plate temperature to be coupled to that of the fluid only through the thermal 
storage term of the plate. And in the steady state, the thermal storage capacity of the 
plate is immaterial since the plate temperature does not change with time. Side condi- 
tion (5) shows this also since in the steady state 90/3t = 0, and the solution must be 



Figure 4. - Influence of coupling parameter « on heat-transfer response to step change in inlet tem- 
perature. 
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independent of r and therefore of e. Also the effect of increasing e on the nondimen- 
sional flux is decreased with increasing rj, as would be expected from the prior 
discussion of the limiting case r\ — 00 . 

A physical interpretation of the coupling parameter e can be made as follows: 
From its definition and the definition of r, 


e 


r “f u 


After some algebraic manipulation this becomes 


e = 



2p C b 
s p, s 


(29) 


But for steady, laminar, slug flow over an isothermal plate, the thermal boundary layer 
thickness is proportional to the square root appearing in the numerator of equation (29). 
That is, 



Therefore, 


e „ p f C P,f 5 t 


Thus, e is a measure of the ratio of the thermal storage capacity of the boundary- layer 
fluid per unit length to the thermal storage capacity of the solid plate per unit length. 
Hence, when e is very small, e « 0, the plate temperature excess will remain at 6 C , 
giving the previously discussed limiting case of Q w = 1. At the other extreme, e ap- 
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proaches «° , meaning that the plate has essentially zero thermal storage capacity com- 
pared to the fluid and thus reacts to impressed thermal loadings instantly and without 
hindering the fluid. This case causes a rather degenerate transient in which is 
zero for t < 1 and equal to its steady- state value for all r > 1, the steady- state value 
being the function of ri given by allowing e to approach °° in equation (17). Doing this 
gives 


% - V” i e ” erfc[Tj] as e — 00 for r > 1 

For this case, once r is greater than unity, the deviation of from unity results 
solely from the thermal history effect. (Recall that the definition of requires that 
it be unity for an isothermal surface at temperature excess 6 . ) An interesting obser- 
vation from figures 2 and 4 is that the nondimensional flux never achieves the infinite or 
extremely large values which are characteristic of solutions in which there is an abrupt 
step in the wall temperature caused by the pure conduction process ahead of the fluid 
front located at x = 0 at t = 0. This observation agrees with the results presented in 
reference 12, where a transient stagnation- point problem was worked for a step in free- 
stream temperature which required a nonzero transit time to be felt at the solid surface, 
When trying to establish conclusions based on figures 2 or 4, we must keep in mind 
that the surface flux has been normalized by dividing it by the steady- state surface flux 
for an isothermal flat plate at fl c . If we compare the ratio of the instantaneous flux for 
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Figure 5. - Heat-transfer response to step change in inlet temperature when fluid and plate are in an 
initial steady state with their temperature not equal to coolant temperature, y = 0. 50; 

<T 0 - Tod/ (Toe - T c ) - 2.0. 



particular values of tj, e, and r to the eventual steady-state flux at the same rj and e, 
and the picture, of course, is vastly different. For instance, at t? = 0. 50 and t = 1, 
the ratio of the surface flux at r = 1 to the eventual steady-state flux (r — 00 ) is 1. 833. 

In figure 5, one representative curve for the transient flux is given for the situation 
where steady-state conditions originally prevailed. Thus, the surface temperature and 
flux are given by equations ( 20 ) and ( 21 ), respectively, and the inlet temperature is T^, 
when suddenly the inlet temperature is increased to the value Tq. For 17 = 0. 50 and 
(Tq - T 00 )/(T co - T c ) = 2, the ratio of the instantaneous nondimensional flux to the 
flux Q sg which existed in the original steady state is plotted versus the time param- 
eter. As shown in figure 5, gg for all r < 1. Then, when the fluid that was 

at x = 0 when the inlet temperature was changed arrives at any position corresponding 
to r = 1 , the flux increases and then gradually decays to its new higher steady-state 
value . 

The simplest approach to transient convection heat transfer is the quasi- steady 
analysis where we assume that the steady- state relations for the wall heat flux are ap- 
proximately valid at each instant of time if the wall temperature at each instant of time 
is used. For the laminar, constant-property, slug-flow problem, the quasi-steady wall 
flux, q Wj qg is given by 



Equation (30) is derived in appendix D by using Lighthill's approach (ref. 19). 

An energy balance on a length dx of plate using equation (30) yields the following 
equation for the quasi- steady wall temperature excess: 

50 f X , 50 

— +a„ / — -dj = b (0 - 0 ) (31) 

3t J 0 Vx-{ 

where 
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and 



The quasi- steady flux is most easily calculated from equation (31) since the integral 
term is directly related to the flux through equation (30) and 0 w is now known from 
equation (32). The result for the nondimensional quasi-steady flux is 

2 2 2 

qs = e‘( € T +2t ? €T ) + Vjr T]e n {erfjer + tj] - erf[rj]} (33) 

Figure 6 shows, for rj = 0. 5 and e = 1, the quasi-steady flux given by equation (33) 
as the dashed line and the exact transient flux given by equation (17) as a solid line. As 
shown, the quasi-steady approximation is not very good for this choice of parameters 
7] and e. In particular, the quasi- steady solution does not predict the lag time t = 1 
needed before a position on the plate can know that the inlet fluid temperature has 
changed. Lyman, reference 12, notes this also in connection with his transit time. In 



Figure 6. - Comparison of exact and quasi-steady solutions for heat-flux response to step change in inlet temper- 
ature when both fluid and plate were originally at coolant temperature. « = 1; r? > 0.5. 
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Figure 7. - Influence of coupling parameier e on quasi-steady solution for heat-flux response to step change in 
inlet temperature when both fluid and plate were originally at coolant temperature, rj = 0. 5. 


figure 7 the nondimensional wall flux is plotted versus e(T - 1) at r\ = 0. 5. (The part of 
the quasi-steady solution for r < 1 has been deleted. ) Again the solid curve is the 
exact solution (eq. (17)), while the quasi-steady approximation (eq. (33)) is shown as a 
dashed curve for two values of e. As would be expected, the quasi- steady solution ap- 
proaches the exact solution for small values of the coupling parameter e. A compari- 
son between the quasi- steady solutions (eqs. (32) and (33)) and the corresponding exact 
solutions for the same case (eqs. (15) and (17)) yields the conclusion that the quasi- 
steady values for both temperature excess ratio and wall flux will be in error by less 
than 10 percent relative to the exact solutions, over the time range r a 1 if e s 0 . 08 
for rj = 0. 5 and e < 0. 20 for 77 = 1 . There is inherent difficulty in making this com- 
parison for the 7 ] = 0 case since the exact solutions tend to zero with increasing time. 
However, if e < 0.02, the error in the temperature excess ratio and the flux will be 
less than 10 percent for rj = 0 if 0 — g(t - 1 )^ 1 . 8 . 

It is interesting to note that if t in the quasi- steady solution is just replaced by 
t - 1, the quasi-steady solution agrees with the exact solution. Quasi-steady solution 
behavior of this general type was also noted by Lyman (ref. 12) in connection with the 
stagnation- point solution. On this basis it would seem reasonable, as a heuristic at- 
tempt t.o getting a better approximation, to just replace r by r - 1 in quasi- steady 
solutions to the more difficult transient convection problems involving a nonslug velocity 
profile. However, this procedure would not be used if axial conduction in the plate were 
significant, because downstream portions of the plate and the fluid would not have to 
wait until the fluid front that was at x = 0 at t = 0 reached them in order that they 
respond to a change in inlet conditions. 

Finally, a representative calculation of e was made for a turbine vane to see 
whether a simpler quasi- steady approach is valid. The conditions used were T^ = 

1390 K (2500° R), a pressure of 82. 6 N/cm 2 (120 psia), u B = 366 m/sec (1200 ft/sec), 
b = 0. 102 cm (0. 040 in. ), x = 3. 25 cm (1. 28 in. ), and p C « 3. 4 J/(cm 3 )(K) 

p, & 
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(50 Btu/(ft^)(R)). The fluid properties were those of air. This yielded e = 1.94x10“ 5 
Hence, based on the idealized model analyzed herein, it would appear as if a quasi- 
steady analysis would be sufficient for gas turbine blades and vanes. Lyman's compu- 
tations for the stagnation- point region, in reference 12, also seems to point to this con- 
clusion. Further investigation, however, is needed to see how the inclusion of the 
effects of axial conduction and nonslug velocity profiles would modify this conclusion. 


SUMMARY OF RESULTS 

By using two successive Laplace transformations, exact analytical solutions for the 
transient gas temperature distribution and surface heat flux are found for the case where 
the gas, flowing over a cooled vane or blade, experiences a step change in its inlet tem- 
perature which then causes both the gas temperature and the vane or blade temperature 
to vary with time. 

1. Equations and curves are presented from which the temperature distribution 
within the gas and the surface heat flux can be calculated. 

2. Quasi- steady results are also derived and compared with the exact transient 
solution, yielding a criterion for use of quasi- steady results with less than 10 percent 
error. 

3. The theoretical model predicts, in this type of transient, the absence of the in- 
finite or very large surface heat fluxes associated with transients initiated by a step 
change in the wall temperature . 

4. A criterion is derived which gives the time needed for the wall flux and the wall 
temperature to- reach a value within 5 percent of their steady- state values. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 14, 1972, 

501-24. 
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APPENDIX A 


SYMBOLS 

a* value which a time parameter, e(r - 1), must be greater than for solution to be 
within 5 percent of steady state 



b plate thickness 

b„ h „/p„CL b 

o c ' ^s p, s 

C constant- pressure specific heat 

r 

c (h c /rk f ) + (yu^/^rajjx) 



coolant- side surface coefficient of heat transfer 

thermal conductivity of fluid flowing over upper surface of plate 

Laplace transform parameter 

instantaneous nondimensional surface heat flux 

same as but for transient initiated after steady state has been reached 

eventual steady- state value of 

quasi- steady equivalent of 

instantaneous surface heat flux 

quasi-steady instantaneous surface heat flux 

P S C P, s b Af 

Laplace transform parameter 
temperature 
coolant temperature 
an inlet temperature 
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an inlet temperature 
t time 

u(t - 1) unit step function, equals 0 for r < 1, equals +1 for r > 1 
u^ free- stream velocity 

x space coordinate along plate 

Y nondimensional y- coordinate 

y space coordinate perpendicular to plate 

o>£ thermal diffusivity of fluid flowing over plate 

/3 expression defined by eq. (B12) 

y temperature excess defined by eq. (B7) 

6^ thermal boundary- layer thickness 




temperature excess, T - 
coolant temperature excess 
steady- state value of 9 
value of 9 at wall 

temperature excess defined by eq. (24) 

T o- 

density 

dummy variable for time 
dummy variable for x 
(h c / k f ) + rp 


V (h c / k f) V «f x / u oo 

t nondimensional time, u^t/x 

Subscripts : 
c coolant 

f properties of fluid flowing over plate 
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qs quasi- steady conditions 
s plate properties 

ss steady state 

w wall conditions 

w, ss wall conditions at eventual steady state 
<*> free- stream conditions 

Superscripts: 

“ function in first transformed plane 

= function in second transformed plane 



APPENDIX B 


DETAILS OF THE INVERSION OF EQUATION (13) OF THE ANALYSIS SECTION 


Equation (13) gives the temperature excess distribution in the second transformed 
plane as 


e =. 


Ue ° ec \iT + rp e 


-yV^ooS+p)/^ 


P(UK + p) 


Ko s + P h 


e c 

S ( U CO S + P) 


+ — + rp 


Qf* 


After defining 


a = — + rp 


(13) 


(Bl) 


and rearranging, equation (13) becomes 



Now the inverse Laplace transformation of equation (B2) with respect to x must be 
found. The second term is found to be 


1 

-pte/Uoo ) 

1 - e 

;" s ( s+ t) 

P 


(B3) 


By use of a table of transforms and the substitution theorem, we get for the first term 
of equation (B2) the following: 
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Thus, equations (B4) and (B3) combine to give the temperature distribution function in 
the first transformed plane as 



Now equation (B5) must be inverted so that p is mapped back into time t. The in- 
verse of the third term in equation (B5) is 



-p(x/u 00 ) 

1 - e 

P 


= 1 - u(t 



where u{t - (x/u^ j] is the unit step function. 

The inverse of the first term of equation (B5) yields the following: 



(B6) 


(B7) 


In attempting to invert the second term of equation (B5), it is advantageous to add and 
2 2 

subtract y u w /4 a^x to the exponent oy+(a^xa J/u^ in order to get e to a power which 
is the square of the argument of the complimentary error function it multiplies, so as to 
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easily perform the inverse transform . Doing this and then defining for convenience 
yield 



rk^ 2r a^x. 


and 



The second term of equation (B5) becomes 


y j e ' p(x/u “ > e d 2 (p«) erfc (- d(p + c) j 


This result can be inverted by using the convolution theorem . For the first factor we 
have 



-y 2 u oo/ 4 a f x 


P 



U oo/ 4Q! f X 


(B8) 


The inverse of the second factor is found by simultaneously employing the substitution 
and translation properties of the transform as 



(B9) 

By using equations (B9) and (B8) and the convolution theorem, the inverse of the second 
term of equation (B5) becomes, using A as a dummy variable for time and denoting the 
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entire second term as n, the following 


if" 1 [II] = - 

P-t 


-Uco y /4Q! f x / 

3 U|t - — 


>fn d 



r- 


fX-iLV / 4d 2 
u 


+ c[X-(x/uj] 


dX 


(BIO) 


The integration is most easily performed if we complete the square on the exponent in 
the integrand by adding and subtracting the term 

h c “f* V u ® y2 
k 2 u °o k f 4a f x 

Equation (BIO) becomes, when c and d are replaced by their equivalents, 


! [( h c o, f x / k f u ») + < h cy/ k f>] / _ 


ifV] =- 

p— t 


u_ 




U, 


/ 

•'x/u. 


1 U co 

exp | - - 

|2r \ a jX 


X--2L + 


2 OjXr /h c yu o 


02 


U oo U oo \ rk f 2r0t f*J 


dX 


(Bll) 


This equation suggests the following variable change: 


J8 = — 

2r 




2a.xr 

JL + J— 

U oo U oo 


X yu 0 


trkj 2ra^xy 


(B12) 


1 / u «> J 
d;3 = — + dX 

2 r V QijX 


(B13) 
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When the quantities 77 , e, t, and Y defined in the analysis are used, the changed 
limits on the integral are as follows: 

X = x/u^ corresponds to 


0 = V + Y 


X = t corresponds to 


0 = e(T - 1) + 77 + Y 


When equations (B12) to (B15) are used, equation (Bll) becomes 


2 / \ /»€(t-1)+t?+Y o 

- -e^W U(T - 1) (A)jj .-r 


e‘ p d/3 


But this, by the definition of the error function, reduces to 


(B14) 


(B15) 


(B16) 


2 

jsr“ 1 [n] - -e 27?Y+7J u(t - l){erf[e(r - 1) + 77 + Y] -erf [77 + Y]} (B17) 

P~t 


After changing to the variables 77 , Y, e, and r in equations (B7) and (B6) and adding 
equations (B6), (B7), and (B17), we arrive at the nondimensional temperature excess 
distribution in the physical plane given by equation (14) of the analysis. 
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APPENDIX C 


SOME DETAILS OF SOLUTION FOR CASE WHERE INLET TEMPERATURE IS 
STEPPED AFTER A STEADY STATE HAS BEEN REACHED 


When fluid at temperature is flowing over the plate of figure 1 (coolant tem- 
perature being T ) in the steady state and then the inlet temperature is abruptly changed 
to Tq, the mathematical statement of the problem is given by equations (1), (22), (23), 
(4), and (5). 

A solution of the following form (eq. (24)) is sought: 

0(x,y,t) = ^(x,y,t) + 0 gs (x,y) (24) 

where is the function given by equation (19) and satisfies the following conditions: 

x = 0 y > 0 6 SS =0 

y — 00 x > 0 0gg “*■ 0 

d6 ss h c , 

y = o x > o -**=-£.(9 -9 ) 

dy k f ss c 

Substituting equation (24) into (1) (noting that 9 co satisfies the steady-state version of 

bb 

(1)) gives the partial differential equation for <p as 



!£ + u = 

9t °° 0x f 3y 2 


(C2) 


The side conditions on (p become 

t = 0' x > 0 y > 0 cp = 0 (C3) 

x = 0 t > 0 y > 0 ^ =T o" T s = ^0 

y — °° t>0 x>0 <pis finite (C5) 
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y = 0 t > 0 x > 0 


(C6) 


+ T^P- 
dy kj 0t 

Now the following transformation is introduced in order to force the problem posed by 
equations (C2) to (C6) to look like the problem given by equations (1) to (5), whose solu- 
tion has already been determined: 

y(x,y,t) = (p(x,y,t) - (p Q (C7) 

Inserting this into equation (C2) gives 


-^ + u 

at 


oo 


3x 3y 2 


(C8) 


The side conditions on y become 


t = 0 x > 0 

y > 0 y = -<p o 

(C9) 

x = 0 t > 0 

y > 0 y = 0 

(CIO) 

y - 00 t > 0 

x > 0 y is finite 

(C11) 

y = 0 t>0 x > 0 

11 — p 0>l + r - 

a y k f L U J at 

(C12) 


Inspection of equations (C8) to (C12) shows that they are identical with equations (1) to 
(5) if 8 q is replaced by -<p q. Hence, the solution function for y is given by equa- 
tion (14) with 0 Q replaced by -(p q, and the solution to the total problem of interest be- 
comes equation (26). 
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APPENDIX D 


DERIVATION OF QUASI-STEADY WALL FLUX EXPRESSION FOR ARBITRARILY 
VARYING WALL TEMPERATURE UNDER SLUG-FLOW CONDITIONS 

The governing equation and side conditions, in terms of 0 = T - T^, for this quasi- 
steady problem are as follows: 


u 


oo 


3x 3y 2 


x = 0 y > 0 0=0 

y — °o x>0 6 is finite 


(Dl) 


(D2) 

(D3) 


y = 0 x > 0 '0 = 0 w (x) (D4) 

Using the Laplace transformation on equation (Dl) with respect to x and also on 
side conditions (D3) and (D4), solving the resulting equation, and applying the boundary 
conditions in the transformed plane yield the transformed temperature excess distribu- 
tion as 


0 = 0 e 
w 


-y V^oP/^f 


(D5) 


where p is the transform variable. Following Lighthill (ref. 19) only the surface heat 
flux in the physical plane will be found by inverting the transformed surface heat flux. 
The transformed quasi-steady surface heat flux is, after some rearrangement for ease 
of inversion, 



(D6) 


Inverting equation (D6) with the aid of the convolution theorem gives the following general 
expression for the quasi- steady surface heat flux for wall temperature varying arbi- 
trarily with x and time t: 
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where i is a dummy variable for x and the integral must be interpreted in the 
Stieltjes sense. 


32 



REFERENCES 


1. Sparrow, E. M. : Unsteady Stagnation- Point Heat Transfer. NASA TN D-77, 1959. 

2. Siegel, R. ; and Sparrow, E. M. : Transient Heat Transfer for Laminar Forced 

Convection in the Thermal Entrance Region of Flat Ducts. J. Heat Transfer, 
vol. 81, no. 1, Feb. 1959, pp. 29-36. 

3. Siegel, Robert: Transient Heat Transfer for Laminar Slug Flow in Ducts. J. Appl. 

Mech. , vol. 26, no. 1, Mar. 1959, pp. 140-142. 

4. Cess, R. D. : Heat Transfer to Laminar Flow Across a Flat Plate with a Nonsteady 

Surface Temperature. J. Heat Transfer, vol. 83, no. 3, Aug. 1961, pp. 274-280. 

5. Goodman, T. R. : Effect of Arbitrary Nonsteady Wall Temperature on Incompres- 

sible Heat Transfer. J. Heat Transfer, vol. 84, no. 4, Nov. 1962, pp. 347-352. 

6. Chao, B. T . ; and Jeng, D. R. : Unsteady Stagnation Point Heat Transfer. J. Heat 

Transfer, vol. 87, no. 2, May 1965, pp. 221-230. 

7. Sparrow, E. M. ; and Siegel, R. : Unsteady Turbulent Heat Transfer in Tubes. J. 

Heat Transfer, vol. 82, no. 3, Aug. 1960, pp. 170-180. 

8. Chambre, P. L. : Theoretical Analysis of the Transient Heat Transfer into a Fluid 

Flowing over a Flat Plate Containing Internal Heat Sources. Heat Transfer, 
Thermodynamics, and Education. Harold A. Johnson, ed. , McGraw-Hill Book 
Co., Inc., 1964, pp. 59-69. 

9. Perelman, T. L. : On Conjugated Problems of Heat Transfer. Int. J. Heat Mass 

Transfer, vol. 3, no. 4, 1961, pp. 293-303. 

10. Soliman, Moustafa; and Chambre, Paul L. : On the Time- Dependent Leveque Prob- 

lem. Int. J. Heat Mass Transfer, vol. 10, Feb. 1967, pp. 169-180. 

11. Inouye, Kenji; and Yoshinaga, Takashi: Transient Heat Transfer near a Stagnation- 

Point. J. Phys. Soc. Japan, vol. 27, no. 3, Sept. 1969, pp. 771-776. 

12. Lyman, F. A. : Heat Transfer at a Stagnation Point When the Free-Stream Temper- 

ature is Suddenly Increased. Appl. Sci. Res. , Sec. A, vol. 13, 1964, pp. 65-80. 

13. Sell, M. G. , Jr. ; and Hudson, J. L. : The Effect of Wall Conduction on Heat Trans- 

fer to a Slug Flow. Int. J. Heat Mass Transfer, vol. 9, Jan. 1966, pp. 11-16. 

14. Rotem, Zeev: The Effect of Thermal Conduction of the Wall Upon Convection from 

a Surface in a Laminar Boundary Layer. Int. J, Heat Mass Transfer, vol. 10, 
Apr. 1967, pp. 461-466. 


33 



15. Davis, E. James; and Gill, William N. : The Effects of Axial Conduction in the Wall 

on Heat Transfer with Laminar Flow. Int. J. Heat Mass Transfer, vol. 13, Mar. 

1970, pp. 459-470. 

16. Luikov, A. V. ; Aleksashenko, V. A. ; and Aleksashenko, A. A. : Analytical Meth- 

ods of Solution of Conjugated Problems in Convective Heat Transfer. Int. J. Heat 
Mass Transfer, vol. 14, Aug. 1971, pp. 1047-1056. 

17. Viskanta, R. ; and Abrams, M. : Thermal Interaction of Two Streams in Boundary- 

Layer Flow Separated by a Plate. Int. J. Heat Mass Transfer, vol. 14, Sept. 

1971, pp. 1311-1321. 

18. Kays, W. M. : Convective Heat and Mass Transfer. McGraw-Hill Book Co. , Inc. , 

1966. 

19. Lighthill, M. J. : Contributions to the Theory of Heat Transfer Through a Laminar 

Boundary Layer. Proc. Roy, Soc. (London), Ser. A, vol. 202, no. 1070, Aug. 7, 
1950, pp. 359-377. 

20. Churchill, Ruel V. : Operational Mathematics. Seconded., McGraw-Hill Book Co. , 

Inc. , 1958. 

21. Abram owitz, Milton; and Stegun, Irene A. , eds. : Handbook of Mathematical Func- 

tions with Formulas, Graphs, and Mathematical Tables. Appl. Math. Ser. 55, 
National Bureau of Standards, June 1964. 

22. Riley, N. : Unsteady Heat Transfer for Flow over a Flat Plate. J. Fluid Mech. , 

vol. 17, pt. 1, Sept. 1964, pp. 97-104. 


34 


NASA-Langley, 1972 33 E-7090 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON. D.C. 20546 


OFFICIAL BUSINESS 

penalty for private use $300 SPECIAL FOURTH-CLASS RATE 

BOOK 


POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 
451 



POSTMASTER : 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 


"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving iimited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. Also includes conference 
proceedings with either limited or unlimited 
distribution. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include final reports of major 
projects, monographs, data compilations, 
handbooks, sourcebooks, and special 
bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from: 


SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 
NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


Washington, D.C. 20546 




